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Abstract. This is an expository note on useful expressions for the density 
function of a product of independent random variables where each variable 
has a Beta distribution. 



1. Introduction 



< 

u 

"t^ . This is a brief exposition of some techniques to construct density functions 

"* (u ■) 
with moment sequences of the form Y[ 7 — r^i where (a)„ denotes the Pochhammer 

symbol ^ia) • ^^^h a density / (x) can be expressed as a certain Meijer G- function, 

that is, a sum of generalized hypergeometric series, and as a power series in (1 — a:) 

whose coefficients can be calculated by a recurrence. The former expression is 

r^ I pertinent for numerical computations for x near zero, while the latter is useful for 

X near 1. 

All the random variables considered here take values in [0, 1], density functions 
JnI . are determined by their moments: for a random variable X we have P [X < a] — 

-^ , J" / (x) dx ioT < a < 1, and the expected value E (X") — L a;"/ (x) dx is the 

nth moment. The basic building block is the Beta distribution (a, /3 > 0) 

1 „,_i ,. .,3-1 






(1.1) h{a,p-x) = -——x^''{l-x 



H. where i3(a,/3):=n£iEl^, then 



^ fa) 

x"h (a, I3; x) dx = , '" ,n = 0, 1, 2, . . . ., 

thus <^ j^ > is a moment sequence if < u < u (with a = u, /3 = v ~ u). The 

moments of the product of independent random variables are the products of the 
respective moments, that is, suppose the densities of (independent) X and Y are 
/, g respectively and define 



(1.2) /*5(x) = 



(.)^//«)»(f)f, 
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then / * g is a density, P [XY < a] — Jq f * g (x) dx for < a < 1 and 

x"{f*g{x))dx= f x^f{x)dx [ y"g{y)dy. 
Ja Jo 

These are the main results: suppose the parameters ui, . . . ,Um and vi, . . . ,Vm 
satisiy Vi > Ui > for each i, then there is a unique density function / with the 

moment sequence I I -. — r— ; 

(1) if also Ui — Uj ^ Z for each i ^ j then for < a; < 1 



J^ r (uk) y ^ r (w, - u,) ^.^||_^^ r {vj - u,) 



\k 

oo m 



n— A;— 1 ^ ^ '^'■ 

(2) for S := X^^^i (t^i — u-i) there is an (ttz + l)-term recurrence for the coeffi- 
cients {cn} such that 

(1.4) fi-) = f^^f[^Ai-xr'{i+±cAi-xr],o<x<i. 

The use of the inverse Mellin transform to derive the series expansion in (|1.3p is 
sketched in Section 2. The differential equation initial value problem for the density 
is described in Section 3, and the recurrence for ()1.4p is derived in Section 4. 

The examples in Section 5 include the relatively straightforward situation m — 
2 and the density of the determinant of a random 4x4 positive-definite matrix of 
trace one, where m = 3. 

2. The inverse Mellin transform 

The Mellin transform of the density / is defined by 

M/(p)= / xP-'f{x)dx. 
Ja 

This is an analytic function in {p : Rep > 0} and agrees with the meromorphic 
function 



T- 



^^T{u,)T{v,+p-l) 



at p — 1,2,3,... thus the two functions coincide in the half-plane by Carlson's 
theorem. The inverse Mellin transform is 

f[x) = — Mf ip) x-Pdp, 

^'^'^ Jcr—ioo 

for a > 0; it turns out the integral can be evaluated by residues (it is of Mellin- 
Barnes type). For each j and each n = 0,1,2,... there is a pole of Mf{p) at 
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p = 1 — n — Wj ; the hypothesis Ui — Uj ^ Z for each i ^ j imphes that each pole is 
simple. The residue at p = 1 — n — Uk equals 

-■— n r(„,r;:!:^„.-„, nrK-^.-..) 

X lim {p — I + n + Uk)T {uk + p — I) ■ 

p— >1— n — Ufc 



To simplify this we use 



r(a^n)=/("\ =(-1)" "^^"^ 



lim {p + po)T{po+p-n) ^ lim {p + po) 



(a - n)„ (1 - a)„ ' 

r{p + po + l) 



_ r(i) _(^i)" 

(-n)„ n! 

Thus 
(2.1) 

f rr^ = (y\ EM.\ V "^""'^ n ^ ("j " "^-^ V TT (i + "fc-^0» n 

(note (1 + Ufc — Ufe)„ = n!);in fact this is a Meijer G-function (see [4l 16.17.2]). 

3. The differential equation 

The equation is of Mellin-Barnes type: let dx '■= ^i-D '■— xdx and define the 
differential operator 

m m 

T{u,v) = -x\{{D + 2- vj) + \l{D + l-Ui). 

The highest order term is (1 — x) a;'" 9™ and the equation has regular singular points 
at a; = and x = 1. We find 

T {u, v) x" Y^ c„a;" = x'^ ^ c„ <^ - J| (n + c + 2 - -yj) a;"+i + [| (n + c + 1 - -u,) x" 

n=0 n=0 [ j"=l j = l 

oo I m m 

= x"" ^ x" < c„ JJ (n + c + 1 - Ui) - c„_i JJ (?i + c + 1 - Vj] 

n=l y j=l j=l 

m 
+ x''coY\_{c+ 1-Mj). 
J = l 

The solutions of the indicial equation are c = Ui — 1,1 < i < m. Assume Ui — Uj ^ Z 
for i j^ j. Let c = ui — 1 then obtain a solution of T (u, u) / (x) = by solving the 
recurrence 

_{\ {ui- Vj + n) nj"li (m ^ -^j + l)n 

'-'" — \\ r I \''n-i ITT'" / I 1^ ''O- 
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Thus the solutions of T (u, w) / = are hnear combinations of 

Ui - Vi + 1,. . .,Ui - Vm + I 



\Ui - U2 + I,. . .,Ui -Ujn + l 
00 m / I 1 \ 

'^-^ -'■-'- [Ui — Ui + 1) 

ioi 1 < i < m (note the factor (uj ~ Uj + 1)^^ = rt!). 

Lemma 1. Suppose g is differ entiahle on (0, 1], g^^> (1) = for < j < k and 
h (x) := {D + s) g {x) then h^^'> (1) = /or < j < fc - 1. Furthermore if k > 
then for n > 

r-i ri 

x^h (x) dx — {s~ n— 1) I x"g {x) dx. 



Proof. By induction d^^D = xdl+^ + jdl for j > 0. Hence h'-^'> (1) = 
gU+i) (1) + (j + s) g(j) (1). Next 

x"h{x)dx = s x"-g{x)dx+ x"+^g'{x)dx 
Jo Jo 

= s x"'g (x) dx + g (1) — {n + 1) / x"'g{x)dx, 
Jo Jo 

and g (1) = by hypothesis. D 

This is the fundamental initial value system: 

(3.1) T{u,v)f{x)^0, 

/(^■) (1) = 0,0< j < m-1 

Proposition 1. Suppose f is a solution defined on (0,1] of L3. 1\) then for 
n>0 



„1 m , s „i 

/ x^f(x)dx^l[p^ / fix)dx 

Jo ,_-, l^Un Jo 



Proof. For < j < m let hj = HLi (-D + 1 - Ui) /, thus hj+i = {D + 1 - Uj+i) hj 
and by the Lemma hY' (1) = for < fc < to — 1 — j. Also /„ x"hj+i {x) dx — 
— (n + Mj+i) /q a;"/ij (x) dx for < j < m — 1. By induction Jj, x"hm {x) dx — 
{-lrmUir^ + ^^)Io^"fi^)dx. 

Similarly /^ x"x JlLi (^ + 2 - v,) f (x) dx - (-1)'" 117=1 (« + «») /o 2;"+V (^) ^^^^ 
Thus the integral — Jq x"T (u, v) f (x) dx implies the recurrence 



/■I ni 

/ x"+v(x)dx=n 
•^0 i=i 



("^ + "^ ^'x"/(x)dx. 



{vt + n) Jo 
Induction completes the proof. D 

Observe that the coefficients {7^} of the solution X]i=i 7j/i i^) of ^^e system 
are not explicit here, but they are found in the inverse Mellin transform expression. 
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4. The behavior near a; ~ 1 and the recurrence 

First we establish, the form of the density / (x) in terms of powers of (1 — x). 
Lemma 2. For a, /3, 7 > and < x < 1 

Proof. Change the variable of integration t = 1 — s + sx then the integral 
becomes 

(1 _ xf+^-^ I (1 - s (1 - x)T-^ s^-1 (1 - sf-^ ds 

^ ^ r(/3 + 7) '''\ /3 + 7 ' )- 

This is a standard formula, see [3l (9.1.4), p. 239] and is valid in < a; < 1 (where 

|l-x|<l). D 

Set 5 := YaLi ("* -■"*)■ 

Proposition 2. There exists a sequence {c„} such that 

Proof. Argue by induction. For m = 1 we have (see (HHJ) 

r(ui)r(i;i -ui) 



'^^^^ (1-) 



T{ui)T{vi - ui) 



"^-"^-ijl + f;^!— j^(l-.x)4. 



Assume the statement is proven for some to > 1, then g = J^h {um+ii Wm+i — Uj^j^i] •) 
has the moments Y\ ."\" . The convolution integral (see (|1.2p ') is a sum of terms 

i— 1 ^ 

c,./ji-o'"-(f) (1-7) ? 

_ /^ ..^«,„+l-l n _ ^\5+n+«„+i-«„+i-l r ((5 + n) r (l^m+l - Um+l) 

r(d + n + t;„+i -u,„+i) 

/ (5 + n, w™+i - 1 
X 2i^i . , ; 1 - a; 

by Lemma[2] ; and x^^+i^^ = 1 + X^^i 1 — ~ (1 ^ ^)"' ■ Thus the lowest power 

of (1 — x) appearing in g is S + w„i+i — Um+i — 1 which occurs for n = 0. By the 



inductive hypothesis 






r ((5) f^^ T (ui) r (u™+i) r {vm+i - um+i) 
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and SO the coefficient of (1 — a;) ^">+i^"™+i^ in g is 

„ ^(J)^(^;„+l -Ujn+i) _ 1 "t-t r (vj) ^ 

T{6 + v„i+i - u„i+i) r{d + v,n+i - Mm+i) r^ r {Ui) ' 

tliis completes tlie induction. D 

For tlie next step we need to express T{u,v) in tlic form x"^ {1 — x) d^ 

3 = 

ables {zi, . . . , Zm\ given by tfie generating function 



^ - g x^ {ttj — bjx)d^. Recall the elementary symmetric polynomials in the vari 



m 

m—n 



SO eo (z) = 1, ei (z) = zi + Z2 + • • • + Zm, 62 (z) = Z1Z2 + . . . + Zm-iZm and e„ (z) = 
Z1Z2 . . .Zm- Thus 

Y[{D + 1- uO - ^ e„_, (1 - u) {xd^y = Y. «J^'^^- 

Let (xd^f = 'El=o^k,^x'^l then xd, (xd^Y = El=o^k,^ {ix'dl + x'+^d',+^), so 
^fc+i,i = ^fe,j-i + iAk,i- This recurrence has the boundary values ^0,0 == 1; ^1,0 = 
Q,Ai^i — 1. The solution consists of the Stirling numbers of the second kind, 
denoted S {k,i) (see [4, 26.8.22]). Thus 

J2 "j^^^^ ^ J2 ^™-j (1 - ") Zl '^ (■?'' *) ^'^^ 

j"=0 j=0 i=0 

rrt m 

^^x^di^S (i, j) e™_, (1 - u) , 
aj ^^S {i, j) Cjn-i (1 - m) , < j < m. 

In particular flm = 1, Um-i — (™)+'5i (1 ^ ""): and ao = e^, (1 — ?i) = H^i (1 ~ %)• 
Similarly 

771 771 771 

E ^^-^'^^ = n (^ + 2 - ^.o = E e™-j- (2 - ^) (2^5,)^' 

i=o j=i i=o 

771 rTi 

m 

6j = E '5' (i, j) em-3 (2 - w) , < j < m. 

The differential equation leads to deriving recurrence relations for the coeffi- 
cients {c„}. Convert the differential operator T {u,v) to the coordinate t = 1 — x; 
set dt ■— ^ (so that dt = —dx)- Write (expanding x^ = (1 — ty with the binomial 
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theorem) 

m m 

T {u, v) = -xY^ bjx^di + ^ a.x^di 
j=o j=0 

rn j + 1 / • I -1 \ m_ J^ / ■\ 

- - E E r t ) ^-'^ '^ (-^*)' + E E ( ■ ) (-ty a, i~o,y 

fe=-m j=max(0,-fe) *- ^-^ ^ ^''^ ^ ^ 

The highest order term (k — 1) is J2T=obj^''^^^i- The term with fc = is 
^ ■ Q (flj — (j/ + 1) bj) P df . The two bottom terms (k = —m, 1 — m) are 

(-ir(a„-6„)ar = o, 

(-1)"-^ (a™_i - &„_i - tdt) dr-' - (-1)"-^ [f2 («4 - zi4) - m - ta, ) a, 

the remaining terms (— m < —k < 0) are 



m— 1 . 

t I 

4 = 1 



Apply T (u, w) to n (with the aim of finding a solution t'' Xl^o '^"^" **-" ^ ("' ''^) / 

m—1 



0); note diV ^ (-f)^ (-7) - ^-^ then 



m 

^-i(7)-E^^(-l)-''(-7).' 

Rrn-l (7) = {-l),n-l [J2{v^-U,)-^-l]= (-7)™-! ('^ - 1 - 7) , 

and 

Rk{i) = {-i)kR'kii), 

These sums can be considerably simplified (and the Stirling numbers are not needed) . 
Introduce the difference operator 

Vg (c) = 5 (c) - 5 (c - f ) . 
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This has a convenient action, for j > and arbitrary k 

W{k- c)^ ^{k- c)j - {k + 1- c)j= {k- c- {k + ] - c:)}{k + l- c)^._, 
= -J (fc + 1 - c)^._i , 

Define the polynomials 



m 

pi ' 



' (c) = n (c + 1 - M,) , g (c) = [] (c + 2 - i;,) 

4=1 i=l 

gi (c) = (1 + cV) g (c) = (1 + c) g (c) - eg (c - 1) . 
Proposition 3. i?_i (7) = gr (7) and for < fc < m - 1 

i?U7)4v'^W7)-(^VV(7). 

Proof. By construction Ejlo Oj*^^**'' = P^^) ^^d Y^^lob^Pdlf = 9(7). 
Apply ■^V'^ to both sides (V acts on the variable 7) of 

m ni 

J=0 i=0 



to obtain 



j—k 

i=fe ^•' ' 



Also 

91 (7) = ( 



m m 

1 + 7V) g (7) - (1 + 7V) J2 f'jt'^tt' - (1 + 7V) J2 (-1)' ^J- (-7), 

j=0 j=0 

= E(-i)'^4("^)^-"^'^(^"^)^-i}^E(-i)'^^-(i+^')(-7).- 

Apply TfeTw^'^ to both sides to obtain 

m .| 

(Fnyi^'" W ^ UTI)! E <-'>'" '^ " + ^■' tFT)! <-■' + «^-' 

i=fc ^-^ ^ 

This completes the proof. D 
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Hence 

oo oo 7n — 1 



-fe 



n=0 n=0 fc=-l 

oo 771—1 



71— 1— m k—-~l 

The recurrence for the coefficients for n > ?7i is 

inin(m,n) 

R„i-1 (n + c) Cn = - X Rm^i-k{n- k + c)cn-k, 
fe=l 

i-^ - c)777-i {5-l-n-c)cn = 

m.in(77i — l,n) 

- X (-71 + fc -c)„^_^_j.i?/,„_;^_j, (n- A; + c)c„_fc -g(n- 771 + c)c„_„, 
fe=i 

where q = for i < 0. At 77 = the equatfon is {—c)^_^ (5 — 1 — c) cq = 0. Let 
c = then for < 77 < 777 — 2 the equations are 

77 

(-n)„_i ((5 - 1 - 77) C„ = - X (-77 + /c)„,_i_;, Kn-l-k (^ " A:) C„_fc, 

fe=l 

but (—77 + k)^_^_i, = for m—l — k>n~k (and < fc < 77) thus the coefficients 
Co, ci, . . . , Cm-2 are arbitrary, providing 777 — 1 hnearly independent solutions to 
T [u, v) f = 0. The recurrence can be rewritten as 



^ m—l 



-1 - 


(777 


-1)!(5- 


-777) 






-1 


(m- 


C71 


(s- 


-1-77) 


\h 



X(-l)'=(7n-l-fc)!i?:„_l_,(77-fc) 



^n — ki 



k=l 



-I f'""^ 1 1 



q{n- m) Cn-m \ ,n>m. 



Assume that (5 ^ Z to avoid poles. But these are different from the desired 
solution which has c = 5 — 1 as was shown in Proposition [51 The recurrence 
behaves better in this case. Indeed 

C77 = - > -, 7— -T i?„_i_fc(71-fc + (5-l)c„_fc 

n ^^ (-"-'5 +1)777-1 

-g (77 - 777 + (5 - 1) c„_m, 



77(-r7-(5+l)„,_^ 



which simplifies to 

inin(77i — 1,77) 

(4-1) c„ = i g (-.-,+ l)/ »--^-^-("-^ + ^-^)^"-^ 
H -, r~~ri 9 {n-7n + 5 -1) c„_„. 

r7(-77-(S+l)„_i 

The term with Cn-m occurs only for n > m. The denominator factors are of the 
form ((5 + 77, — 1) ((5 + 77 — 2) ... ((5 + 77 — fc). If 77 > 777 then the smallest factor is 
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S + n — ni > S > 0; otherwise the smahest factor is 6 (for k = n). Hence this 
sohition is well-defined for any S > 0. 

Theorem 1. Suppose ui, . . . , u,„, ui, . . . , Vm satisfy Vi > ui > for each i then 
there is a density function f on [0, 1] with moment sequence Yl ( "'"'^ 



Z— 1 



/w = r^nfgya-)-{i + i:-a 



=1 



where the coefficients {c„} are obtained with the recurrence (RTll) using cq — 1, and 

S = YT=li'"^-Ui)■ 
PROOF. The density exists because it is the distribution of the random variable 
Yi^i -^i where the X^'s are jointly independent and the moments of Xi are ;"\" 
for each i. By Proposition [2]/ has the series expansion given in the statement. Let 
g (x) be the function given in the statement and suppose for now that 5 > m then g 
is a solution of the differential system (|3.1I) (because of the factor (1 — x) ^ ). By 
Proposition [1] (7 has the same moments as C f for some constant C . By Proposition 
[2]/ and g have the same leading coefficient in their series expansions. Hence f = g. 
The coefficients Cn are analytic in the parameters for the range S > m. Each 
moment /„ a;"/ (x) dx is similarly analytic and so the formula is valid for all (5 > 0, 
by use of analytic continuation from the range 5 > m. D 



The coefficients occurring in the recurrence (j4.ip are expressions in the pa- 
rameters u, V, which can be straightforwardly computed, especially with computer 
symbolic algebra. 



5. Examples 

5.1. Density for m — 2. For the easy case to = 2 we can directly find the 
density function, in a slightly different form. 

Given ui,U2,vi,V2 > and d = vi + V2 — ui — U2 > set 



r(z;i)r(«2) 

g{u,v;x) - 



then 

f-i 



F (ui) F {u2) F (S) 



a:"5 (u, v; x) dx = " " :^ = 0, 1, 2, 

(«l)„ (W2)„ 



PRODUCTS OF BETA VARIABLES 11 



Proof. Consider 



I ^rn-«i i (1 _ a;) 2F1 I ■,l-x\dx 

{v2 - U2)„ {vi - U2)„, r (n + Ml) T{5 + m) 



00 



E 



"^' {^)m r (ui + (5 + m + n) 



r (Wl) (Ul)„ r ((5) y^ (i;2 - U2)„j (wi - U2), 



r(n + ui+J) ^^ m\{n + ui + 5)^ 

_ r (-m) (mi)„ r (^) r (n + Ml + (5) r (n + Ml + (5 - t;2 - t'l + 2U2) 
T {n + ui+ 5) T {n + ui+ 5 — V2+ M2) F (n + ui + 5 - ui + 1*2) 

_ r (mi) r (M2) r (1^2 + ^i - tti - ^2) (mi)^ (•»2)n 

r(i;2)r(z;i) {v2)Avi)J 

for each n. D 



By using standard transformations we can explain the other formulation for 
I {u, v; x) near x = Q. From ^ p. 249, (9.5.7)] 

rl""^^;! x\ r(7-a-^)F(7) ^/ a,^ .^ 

7' / F(7 — q:)F(7 — (5) \l + a + (5 — 7' 

r(Q; + (5-7)F(7) ^_a-5i:. /^ 7-a,7-'5 



^ r(a)F(5) Vl + 7-a-^' 



applied to g (u, w; a:) (provided u\~ U2 ^ Z) wc find 

F(w2)F(wi)r(u2 -ui) 



(5.1) g{u,v\x) 



+ 



F (wi) F (U2) F (w2 - Ml) F (ui - Ml) 

1 - X) 2-Fl 

V 1 +■! 

F(m2)F(mi)F(mi - U2) 



x.--(i-.)^-SFirv""''""^- 

\ 1 +Ui - U2 



F (mi) F (U2) F {v2 - U2) F (wi - U2) 



x.--(i-.)^-SFirv""''""^- 

\ 1 + U2 - Ul 

This is quite similar to the general formula p.ip , and the following standard trans- 
formation explains the difference 

(5.2) 2i^i ' ;a:; =(l-a;)' " ' 2i^i ' ;x 



If Ml — 1*2 G Z then there are terms in logx. The relevant formula can be found in 
[H p. 257, (9.7.5)]. Suppose U2 = ui + ^^ and n = 0, 1, 2, . . . , S — V2 + vi — 2ui — n 
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then 



/ N r(w2)r(wi) s-i 

g (U, V; X) = r r r (1 — X) 

r (ui) r [ui +njT [V2 - ui) r [vi - ui) 

-1 v^ (n — fc— 1)! J. 

^^ ' Z^ fci ("2 - ui - n)fc {vi - wi - n)^ (-a;) 

fe=0 

r (mi) r (ui + 7i) r (w2 — ui — n) r {vi — ui — ?i) 



xi-lrx^-+--'i-\ogx)-2F^^ ^ ^' ;x 

n! \ n + 1 

(-l)"r(z;2)r(t;i) ni+„-l ., _ ^^~l 



r (mi) r (lii + n) r (u2 — ui — n) r (ui — ui — n) 

(l^2 - Ml)fc (I'l - Wl)fe 



(1-^r 



X 



E 



..0 ^!(- + ^)! 



fc 



{V' (fc + 1) + V' (" + ^ + 1) ~ V^ (^2 — wi + ^) — ^ (t^i — ui + A:)} x 



li ui — U2 ^ Z then near a; = the density is ^ Coa;"^ '^ + Cia::"^ ^, but if 
M2 = ui + ?i then the density - Coa:"!-^ + Cix"i+""i (- log a:). 

5.2. Example: parametrized family with m — 3. Consider the determi- 
nant of a random 4x4 state, that is, a random (with the Hilbert-Schmidt metric) 
positive-definite matrix with trace one. The moments can be directly computed 
for the real and complex cases and incorporated into a family of variables with a 
parameter. Here the variable is 256 times the determinant (to make the range [0, 1]) 
and a = h for M, a = 1 for C, and a = 2 for H (the quaternions). This example is 
one of the motivations for the preparation of this exposition. The problem occurred 
in Slater's study of the determinant of a partially transposed state in its role as 
separability criterion 5^. 

The moment sequence is 

(!)„(« + 1)„(2« + 1)„ ^_o,l,2,...; 



{3a+l)j3a+l)j3a+l) 

thus 6 — 6a + ^. For generic a the density is 

3(12a + l)(6a + l)(4a + l) /| - 3a, i - 3a, i - 3a 

64a2 ^ \ 1 - a, 1 - 2a ' ^ 

r (6a + I) r (3a +1)210" ^ ^ - 2a,'^ - 2a,l - 2a ^ ^ 



4asin(7ra)r(a + l)r(8a + l) V 1 - a, 1 + a 

„2a 



(2a -I- 1)^ TT^r (3a + |) T (6a + f ) 2-^" 



48 sin (Tra) sin (2™) T (2a -I- i) T (a + §) T (a -I- 1)" 



3 -a i-a i 

'2 '-^'4 

l + a,l + 2a 



X 



X 3i^2( * -"'2 ,"'f ".^ 
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For numeric computation at a = |;1)2 one can employ interpolation techniques; 
for example 

2 1 

/ (ao; x) = -{f (ao + h;x) + f {ao - h; x)) - g (/ ("o + 2/i; x) + f (ao - 2h; x)) 

where / (a; x) denotes the density for specific a and the last term is the error (for 
some ^ G (—2,2)); thus the perturbed densities can be computed by the general 
formula. 

5.3. Example: the recurrence for m = 4. Given ui, . . . , U4, vi, . . . ,V4 de- 
fine p(c) = nLi(c+l-'«0,g(c) = nLi(c + 2-w»), gi(c) = (c+f)g(c)- 
cg (c - 1), (5 = J2t=i («» - "0 : 

K (7) = P (7) - 91 (7) , 

^'i(7)=Vp(7)-^Vgi(7), 

^2(7) = iv2p(7)-ivV(7), 
then set cq = 1, 

ci = T^2(^- l)co, 

^^-2(m)^^(')^^-^^^'^('-^)^°' 

^^ - 3(^^^ (^ + ^) ^- 3(mW^^'^ (') ^^ 
^ i?^ (,5-1) CO, 



—-J- — i?2 {n + S-2) c„_i — --R[ {n + S-3) c„_2 

n (0 + n — f ) rt (d + n — 2)^ 

R'q (n + (5 - 4) c„_3 H — — — -g (n + (5 - 5) c„_4, 



n(^ + n-3)3 "' ' n{5 + n-3) 

for n > 4. 

5.4. Example: a Macdonald-Mehta-Selberg integral. Let 5" be the 3- 

dimensional unit sphere -^ a: G M"* : X]i=i ^? ~ 1 r with normalized surface measure 

duj. Consider ni<i<7<4 (xi ~ Xj) as a random variable (that is, evaluated at a dcj- 
random point). Interestingly, the maximum value j^ is achieved at the 24 points 
with (permutations of the) coordinates < ±^v9 ± 3\/6 L which is the zero-set of 

the rescaled Hermite polynomial if 4 [V^t) . The Macdonald-Mehta-Selberg integral 
(see [U p. 319]) implies (for k > 0) 

•'Sl<^<J<i ^ T{2 + 6K)r{l + K) 
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For integer values k — n the Gamma functions simplify to Pochhammer symbols; 
then by use of formulas like {l)^„ = 4^" (j), (5) (|) (1)„ the value becomes 

,, _ 1 ilJn y2)n ViJn 



108"(i)„(l)„(i) 



Let fn denote the density function of £> = 108 ni<i<7<4 (^i ~ ^j) (^Y the general 
results the range of I? is [0, 1]). Applying Theorem [T] we find 

, , , \/2,, ,1 f 221,^ , 156697,^ ,2 232223093,^ ,3 

fjj(x) ^ -^(l~x)-^ <1 + 1 - a;) H 1 - xY H 1 -xV + . 

■' ^ ' 37r ^ '^ \ 216^ '^155520^ ' 235146240^ ' 

By formula (fOI 



/d [x) = "fix' 



+ 13X~ 




I2X 



3^2 




where 



71 



72 



73 



3r(!)^r(^)r(ii) 

2%/3 



37r 



37r3 



r - 



^'^ 



r ^ 



It is straightforward to derive a series for the cumulative distribution function 
Fd {x) = /q^ Jd [i) dt. Figures [T] and [5] are graphs of fo and Fq respectively (of 
course there is vertical asymptote for fu). For computations we used terms up to 
the eighth power, with the series in x for < x < 0.55 and the (1 — x) series for 
0.55 < X < 1. For a better view there is a graph of Fd (x) for < x < 0.04 in Figl2] 
and of 1 - Fd (x) for 0.4 < x < 1 in Fig. H 
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Figure 1. Density of D, partial view 




Figure 2 . Cumulative distribution function of D 
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Figure 3. Part of cumulative distribution of D 




Figure 4. Part of complementary cumulative distribution of D 



